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Abstract: We compute the expectation values of circular Wilson loops in large 
representations at strong coupling, in the large- limit of the M = 2 superconfor- 
mal theory with SU{N) gauge group and 2N hypermultiplets. Employing Pestun's 
matrix integral, we focus attention on symmetric and antisymmetric representations 
with ranks of order A^. We find that large rank antisymmetric loops are indepen- 
dent of the coupling at strong 't Hooft coupling while symmetric Wilson loops grow 
exponentially with it. Symmetric loops display a non-analyticity as a function of 
the rank, characterized by the splitting of a single matrix model eigenvalue from the 
continuum, bearing close resemblance to Bose-Einstein condensation in an ideal gas. 
We discuss implications of these for a putative large- string dual. The method of 
calculation we adopt makes explicit the connection to Fermi and Bose gas descrip- 
tions and also suggests a tantalizing connection of the above system to a multichannel 
Kondo model. 
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1 Introduction and summary 

Wilson loops, while being fundamental probes of gauge theory physics, following 
the development of the AdS/CFT correspondence [1], are now also known to play 
a key role as probes of string duals of large- A gauge theories [2-4]. This aspect of 
Wilson loops was sharpened by the beautiful results of [5, 6] for supersymmetric BPS 
Wilson loops in A" = 4 supersymmetric (SUSY) Yang-Mills theory. The upshot of 
these works was the conjecture that circular supersymmetric Wilson loops in N' = A 
SYM are computed exactly, to all orders in the gauge coupling, by the Gaussian 
matrix model. This conjecture was subsequently proven by Pestun [7] by considering 
Af = 2 and A/" = 4 SUSY gauge theories on a four-sphere and computing their 
partition functions using powerful localization techniques. 

One of the theories for which Pestun provided a finite dimensional matrix integral 
representation is the SU (A) M = 2 superconformal theory (SOFT) with Nf = 2N 
flavour hypermultiplets (also referred to as the Ai A/" = 2 SOFT). Recently, Passerini 
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and Zarembo [8] explored properties of this matrix integral in the large- iV limit and 
deduced the behaviour of the circular Wilson loop in the fundamental representation, 
at strong 't Hooft coupling. The focus of the present paper is to use the large-A^ 
matrix integral for the M = 2 SCFT to explore properties of supersymmetric Wilson 
loops in large representations, when the rank of the representation is of order A^. 
Our primary motivation is to compare and contrast the results with corresponding 
quantities in the A/" = 4 theory at strong coupling, and hence draw inferences about 
the possible nature of the large- string dual, if any, of the M = 2 SCFT. 

The Nf = 2N superconformal theory has vanishing beta function and therefore 
an exactly marginal gauge coupling. Finding a string dual to this theory in the large- 
N Veneziano limit [9] at strong 't Hooft coupling, is a long standing problem. Recent 
proposals in this direction include [10-13]. In all cases the proposed backgrounds 
either contain regions of high curvature or are partly non-geometric as in [10, 13]. 
It is therefore interesting to ask if field theory probes such as Wilson loops can shed 
light on features of a putative large- iV string dual. Experience with the A/" = 4 
theory indicates that Wilson loops in generic (large) representations can indeed act 
as effective probes of the dual geometry. In particular Wilson loops in symmetric 
and antisymmetric tensor representations are computed by probe D3 and D5-branes 
in AdS^ X [14-21]. These are known to have world- volumes AdS2 x S"^ C AdS^ 
and AdS2 x {S^ C S^) respectively, and thus are sensitive to different aspects of 
the geometry. 

In [8], a careful evaluation of the circular Wilson loop in the fundamental repre- 
sentation was performed at strong coupling in the M = 2 SCFT, making use of the 
corresponding matrix model. This led to a natural identification of a dual effective 
string tension T^a = ^ In A, where A is the (large) 't Hooft coupling. The curious log- 
arithmic dependence on the 't Hooft coupling translates to a non-exponential growth 
of the Wilson loop with A at large A. The non-exponential dependence had already 
been deduced via scaling arguments applied to the associated matrix model in [22]. 
This is to be contrasted with the circular Wilson loop in A/" = 4 SYM which grows 
exponentially at strong coupling as ~ e^. 

A particularly curious feature of the one-matrix model for the Ai JV = 2 SCFT 
is that the large- A^ eigenvalue distribution has an infinite support at infinite 't Hooft 
coupling, with most of the eigenvalues remaining finite^ . The limiting form of 
the eigenvalue distribution (or the large- A^ "master field") does not describe the 
fundamental Wilson loop since the latter diverges in the limit A — oo. However, in 
this paper we show that the limiting distribution does describe Wilson loops in large 
enough tensor representations. 

Specifically, Wilson loops in the antisymmetric tensor representation with rank 

^This was realized in [8]. In earlier work [22], it was noted, using scaling arguments that the 
eigenvalue density approaches a limiting shape. 
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k of order N can converge to a result that is independent of A. In the large- limit 
with / = fixed, we find that the antisymmetric Wilson loops are determined by 
the endpoint of the eigenvalue distribution only for / -C vhiA/A <^ 1, and beyond 
this range approach a regular limit determined completely by the limiting eigenvalue 
distribution at infinite 't Hooft coupling. This behaviour strongly suggests that the 
internal/compact factors in the large- string dual, probed by the corresponding 
D-branes, remain highly curved or non-geometric as in [10]. 

The situation with Wilson loops in the symmetric tensor representation of rank k 
turns out to be somewhat different. Their expectation values are determined by the 
endpoint of the eigenvalue distribution at strong coupling and essentially track the 
behaviour of the fundamental Wilson loop up to a critical value of / = /c ~ Vln A/A. 
Beyond this point the Wilson loop experiences a non-analyticity characterized by the 
splitting of a single eigenvalue from the rest of the large- distribution. A related 
non-analyticity was observed for the symmetric loop in A/" = 4 SYM, in [20]. In that 
case, the position of the split eigenvalue from the large- A^ distribution was mapped to 
the position of the probe D3-brane in AdS^ x [15] which computes the symmetric 
Wilson loop. 

The approach we use to compute the large rank Wilson loops is identical to that 
of [20]. In this paper we further emphasize the connection of the symmetric and 
anti-symmetric representations to the free Fermi and Bose gas pictures. This in turn 
suggests a potential tantalizing connection to the multichannel Kondo model [23-26]. 

The most straightforward inference we can draw from our results is that the 
string dual to the M = 2 SOFT should have a weakly curved AdS^ part which can 
be probed by any D-branes that compute Wilson loops in symmetric representations. 
The exponential growth of the latter with the 't Hooft coupling indicates that the 
corresponding D-branes which compute them must be semiclassical, bearing some 
resemblance to the situation in the A/" = 4 theory. On the other hand, the behaviour 
of Wilson loops in the antisymmetric representation suggests that the internal or 
compact factor of the geometry must be highly curved. Our results may be viewed 
as predictions for the tensions of corresponding probe D-branes. 

2 The matrix model 

In [7] it was shown that circular, supersymmetric Wilson loops in A/" = 4 and M = 2 
superconformal gauge theories on are computed by matrix integrals as a conse- 
quence of localization of the partition function onto constant configurations. For an 
Af = 2 SUSY theory, the supersymmetric Wilson loop, labelled by some representa- 
tion R of SU{N), is defined as 




(2.1) 
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where $/, (/ = 1, 2) are the two real adjoint scalars of the M = 2 vector multiplet 
and is a constant unit vector in M?. Any specific choice for this vector breaks 
the SU{2)r X U{1)r global R-symmetry of the AT = 2 SCFT to SU{2)r. In the Ai, 
Af = 2 SCFT with SU{N) gauge group and Nf = 2N hypermultiplets formulated 
on the four-sphere, the expectation value of the supersymmetric Wilson loop Wr{C) 
is computed by a corresponding observable in a one-matrix model, 

{Wr{C))x=2 = {TiRe^'"')^^. (2.2) 

The N X N matrix "a" can be viewed as the (constant) VEV for one of the adjoint 
scalars in the vector multiplet. In the diagonal basis for a = diag(ai, 02, . . . oat) with 
^ Oi = 0, the partition function on 5"^ is obtained by integrating out all massive 
adjoint modes and the fundamental fiavour fields at one-loop order, and subsequently 
integrating over the diagonal VEVs 

i^RiC)) = uolr^^n^ [ [He^^''''"'2i-ioop(a)|2i„st(a)pTrH,e2™. (2.3) 



Vo\{SU{N)) 

We have suppressed the explicit form of the integral in terms of the eigenvalues 
{a,}. In particular, in the diagonal gauge the measure factor will give rise to a 
Vandermonde determinant. The first term in the integrand arises from the conformal 
coupling of the scalar fields to the curvature of S^. The perturbative contributions 
are 1-loop exact obtained by integrating out fiuctuations about the VEV, and Zjnst 
is Nekrasov's instanton partition function [27] for the equivariant theory on M^. In 
this paper we will be interested in the large-iV limit in which we expect instanton 
contributions to be exponentially suppressed and a careful consideration of these 
terms confirms this expectation [8]. In terms of the eigenvalues, the contribution 
from 1-loop fiuctuations takes the form. 
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The specific form of the one-loop correction is directly related to the vanishing of the 
beta function for the gauge coupling. When the number of hypermultiplets differs 
from 2iV, a divergence ensues in the one loop fluctuation determinant above, which 
then renormalizes the coefficient in front of the quadratic term in the action i.e. the 
gauge coupling. 

We are only interested in evaluating Wilson loops at large- with A = g'^N flxed. 
Therefore, what is needed is the large- A^ eigenvalue distribution p\{x) which solves 
the saddle point equation 

^-x-K{x) =V dypx{y) (^-^^-K{x-y)^ , xe[-^l,p] (2.5) 

/■^ H'ix') 
/ PA(x)rfa: = l, K{x) ^ -— ^. 
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Recall that the spectral density of the matrix model is defined as the large- 'con- 
tinuum' limit of 



P\{x) = ^^^i^-ai). (2.6) 



An exact solution of the saddle point equation is not (yet) known, but crucial proper- 
ties of px{x) can be inferred from the behaviour of function K{x) which appears both 
as a central force term and in the pairwise interaction of eigenvalues. As observed 
in [8], for small x, K{x) ~ 2({3)x^ while for large x, K[x) — > 2xlnx. This implies 
that the pairwise interaction between eigenvalues is repulsive at short separation 
(dominated by Vandermonde repulsion) and attractive at very large eigenvalue sep- 
aration. On the other hand the central quadratic potential (attractive) dominates 
at short distances, but is overwhelmed by the repulsive K{x) at large distances. 
Importantly, at large distances the one-body term K{x) precisely counteracts the 
two-body force K{x — y), so that for a large enough spread of the eigenvalue distri- 
bution, the behaviour at the endpoints is controlled by the Vandermonde repulsion 
and the quadratic one body potential. 

We list below the main consequences [8] of these observations: 

• The spectral density associated to the eigenvalues of the matrix model (2.3) 
above, attains a limiting form in the limit of infinite 't Hooft coupling, which 
is independent of the coupling. 

The spread of the eigenvalues is infinite (/i — oo). An important feature of the 
limiting form of Poo{x) is that it cannot be used to yield a finite expectation 
value for the Wilson loop in the fundamental representation. In particular, 
^^i^^Ia-^oo ~ I^oc dx poo{x) e^'^^ is divergent at face value. 

• Hence, the finite A corrections to the exponential tail of the limiting distribu- 
tion are crucial for determining the correct value of the Wilson loop at strong 
coupling. For finite (but large) A, since the force on an eigenvalue at large dis- 
tances is determined by the quadratic one-body potential and the Vandermonde 
repulsion, it can be argued that near its endpoints the eigenvalue distribution 
should smoothly interpolate between the limiting distribution and the Wigner 
semi-circle law, 

px{x) ~ ^V/"^ ~ ' a; ~ /i > 1. (2.8) 

The location of the endpoint p can be estimated by requiring the interpolating 
distribution to be correctly normalized i.e., dx p\{x) = 1, assuming that 



-5 - 



the crossover between (2.7) and (2.8) occurs at x ~ This imphes, 



A ~ v^e"^/2 , /i = - hi A + . . . . (2.9) 

TT 

It is then straightforward to infer the A-dependence of (Wci) in the large-A 
regime. The relevant integral is dominated by the endpoint of the spectral 
density, i.e., Eq.(2.8), so that 

(l^n) = i^^. (2.10) 

Comparison with the standard result ^ for the Wilson loop at strong 't Hooft 
coupling, in theories with weakly curved AdS duals [1-6], then suggests that 
the effective string tension is T^s = ^ In A in a putative string dual of the 
Nf = 2N theory. 



3 High rank k ~ 0{N) Wilson loops 

In the context of the AdS/CFT correspondence it has been known for some time 
that Wilson loops in tensor representations of rank k with k ~ 0{N), in the large-A^ 
limit, are computed by semiclassical D-brane probes in the (weakly curved) dual 
geometry [14-18, 20] of the form AdS^ x X^. The antisymmetric and symmetric 
tensor representations are computed by k fundamental strings 'blown up' into probe 
D5- and D3-branes, respectively. The former are probes of the internal geometry 
and their dependence on the ratio is determined essentially by the volume of 
the four-cycle inside X^, wrapped by the puffed up D5-brane. On the other hand, 
the D3-branes are embedded completely in the AdS^ directions with world-volume 
AdS2 X S\ 

It is conceivable that the behaviour of such high rank Wilson loops in the M = 2 
superconformal theory will contain some hints of a large-A^ string dual. It is a priori 
unclear whether such a string dual will have weak curvatures or not, but we expect 
that exact results from field theory may allow us to draw some inferences. 

3.1 Antisymmetric representation 

Given that supersymmetric Wilson loops in the M = 2 SOFT are computed by 
Pestun's matrix model (2.3), Wilson loops transforming in various representations 

^For large-A^ superconformal gauge theories with a supergravity dual (e.g. A/" = 4 SYM) , the 
Wilson loop W ~ ^^.-^{y/XA) / \^/^ . Here A is the area of a minimal string (disk) world-sheet in 
AdS, whose boundary traces out the Wilson loop in the dual gauge theory. The pre-factor of X~^/^ 
arises from gauge-fixing on the world-sheet and depends on the number of zero modes, determined 
by the Euler character of the world-sheet [28] . 
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can be expressed as expectation values of appropriately symmetrized polynomials 
of eigenvalues of the random matrix e^^"". In particular, the rank k anti-symmetric 
tensor representation is explicitly 

^ ^ ^ ^ l<ii<i2<...<ifc<A' 

where the expectation value on the right hand side is within the random matrix 
model (2.3) and the result is normalized by the dimension of the representation, 
dim(A,)=Q. 

To evaluate this, we will follow the technique of [20], assuming that the insertion 
of the Wilson loop operator does not change the large- eigenvalue distribution of 
the matrix model itself. This assumption is consistent in the strict large-A^ limit, 
since the matrix model action is (9(A^^), whilst the insertion of the rank k Wilson 
loop operator introduces a 'perturbation' of order and cannot alter the large-A^ 
saddle point distribution. 

Since our primary interest is in the strong coupling limit of the M = 2 gauge 
theory, we will compute the expectation value (3.1) simply by using the limiting form 
of the eigenvalue density Poc{x)- Despite the fact that this distribution cannot provide 
a finite result for the fundamental representation, we will find that large enough 
antisymmetric tensor representations (-^ fixed as — t- oo) can converge to a smooth 
result (at infinite A), independent of the 't Hooft coupling A. (Independently of this, 
the Wilson loops for such large rank representations will always scale exponentially 
with A^). 

The Wilson loops in the antisymmetric representation can be read off as coeffi- 
cients of the characteristic polynomial for the matrix e^™: 

r df 1 ^ 

(^A.) = / ^^ra(n(^ + ^'™^))— (3-2) 

.7 = 1 



The large- A^ limit provides a particularly convenient representation of this expression 
as an integral over the auxilliary spectral parameter t, which can then be evaluated 
in a saddle point approximation. 



N 1^ dxpxix) ln(l + e"2"^r^) 



(3.3) 



Note that this integral representation takes in the eigenvalue distribution (at 
infinite A^) as input, and can be used, in principle, to evaluate higher rank Wilson 
loops for any A. It also enjoys a symmetry under the operation k ^ N — k, which can 
be understood by performing the variable change t —> and from the symmetry 
of p\{x). This symmetry is equivalent to charge conjugation, so that Wa^. = Wa^ 



^N-k 
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Figure 1: The integral along the contour C, on the cylinder yields the rank k antisym- 
metric Wilson loop. The contour lies to the left of the branch cut between — + | and 
/i + I where // ~ The saddle point on the first sheet lies on the real axis at Zsp, well 

away from the branch cut. 

as expected for sources transforming in the antisymmetric tensor representation of 
SU{N). 

Another feature of the formula (3.3), is that the exponent (as a function of 
t) has a branch cut singularity. This branch cut arises from the zeroes of the 
characteristic polynomial coalescing in the large- iV limit to form a continuum. The 
branch points are at t = —e~'^'^^ ~ — A^''(lnA)^ and t = —e^""^ ~ A^(lnA)"^. In the 
limit of infinitely large 't Hooft coupling the branch cut stretches along the entire 
negative t-axis, and we must be careful whilst considering this limit to evaluate the 
Wilson loop. 

In the limit k, N ^ oo with fixed, the t-integral can be evaluated by the 
method of steepest descent. The condition for the existence of a saddle-point in the 
t-plane is then, 

rdx-^-f = 0, /^A. (3.4) 

An interesting and potentially useful interpretation of the system can be obtained 
by adopting a different parametrization, 

t = . (3.5) 

This maps the plane to the cylinder with a branch cut along — /i + ^ < z < +fi + ^ 
(Figure (1)). Along the branch cut the integrand develops an imaginary part. Given 
the analytic structure of the integrand, we can use this procedure as long as the 
location of the saddle points is away from the branch cut singularity. With the 
exponential parametrization, the saddle point equation acquires a nice interpretation, 

["dx , =/. (3.6) 

This equation gives the occupation number of fermions in the presence of a chemical 
potential z and with a density of states px. Depending on the specific form of px{x), 
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we may be able to infer an effective "temperature" of the system by appropriate 
rescalings of the integration variable. 



3.2 Infinite A limit 

The large-A distribution function has the property that its endpoints run off to 
infinity in the limit of infinite 't Hooft coupling. We know that the fundamental 
Wilson loop (2.10) diverges in this limit because its expectation value is dominated 
by the largest eigenvalues in the vicinity of the endpoint of the eigenvalue density 
Pa- However, we can expect large rank representations to behave differently. While 
individual terms in the sum (3.1) may diverge as A is taken to infinity, in the limit of 
large k, N with fixed, the total number of terms in the sum grows exponentially 
with A^. As the distribution poo is peaked about x = 0, a majority of eigenvalues 
reside in the bulk of the distribution and their contribution to the sum can dominate 
the result for high rank representations. Indeed, this is precisely what occurs. 

The integral representation for the antisymmetric loop at A = oo explicitly reads, 

(^^.)Ia^oo c?^exp[Ar(21n(coshf + ^)-(/-i)27rz + ln4)] . (3.7) 

Despite the fact that the integrand in the original expression (3.3) for Wai, is mani- 
festly periodic under z ^ z + i, its explicit form above in the A — oo limit, is only 
periodic under the shift z ^ z + 4i. This is due to the branch cut shown in Figure(l) 
which now has an infinite extent. 

The large-A^ saddle point equation follows directly from the above result, 

^sinhf (coshf + ^)"' = f-l- (3.8) 

This equation has distinct roots with Imz = and Im^; = 2i. For the moment, we 
focus our attention only on the real root. The complex root lies on a different sheet, 
due to the branch cut discussed above. We find one saddle point on the real axis 
satisfying, 

coth (I ^.p) = (3.9) 

J 2 

As a function of /, in the domain < / < 1, ^sp takes on all values on the real axis. 
For small / and / approaching unity, the saddle point runs of to infinity. 



z, 



sp 



fln/ + ..., /«1, (3.10) 

= -|ln(l-/) + ... , (1-/)«1. 

A nice property of this solution as a function of /, is that 

ZsM) = -Zsp{l-f). (3.11) 
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This reflects the symmetry of the integral representation for the antisymmetric Wil- 
son loop, under z — )■ —z accompanied by the replacement f ^ 1 — f . In the free 
fermion interpretation, ^gp acts as a real chemical potential. Since the distribution 
Poo{x) is independent of A, the effective temperature for the Fermi-Dirac distribution 
in Eq.(3.6) is 0{l). 

We now proceed to evaluate the Wilson loop itself. After deforming the contour 
to pick up the saddle point contribution (formally, before taking the A — ?■ oo limit) 
from the first sheet, we have 



* ' I A— ^oo 



exp 



21n 



sinh (f ^sp) 
f -- 

J 2 



(/-l)27rz 



sp 



(3.12) 



The explicit expression is not particularly illuminating, but the expansion of the 



exponent in a power series in / 



N 



is interesting: 



N 




k_ 

N 



< 1. 



(3.13) 

The most notable features of this formula^ are the non-analytic (logarithmic) 
dependence on / for small /, the complete absence of any dependence on the 't 
Hooft coupling A and that beyond the leading order, subsequent terms are organized 
in a power series in (as opposed to for instance). The logarithmic behaviour 
is an immediate consequence of Eq.(3.12) and that Zgp ~ In / for small /. 



3.3 Large but finite A 

For any finite value of the coupling constant, the eigenvalue distribution has a finite 
extent [8], and the exponential tail of p\{x) smoothly matches on to a Wigner semi- 
circle distribution which vanishes a.t \x\ = p ^ ^In X + . . .. In the fixed k, large- 
limit, we expect that the value for the (exponent of) the antisymmetric Wilson loop 
is k times that of the fundamental loop (from large- A^ factorization). The result 
for the fundamental loop (2.10) depends on A, whereas we have seen that for large 
enough representations with f = fixed, the observables have a A-independent limit 
at large A. We will now see how these two results can be reconciled. 

To understand that there must be a qualitative change in the behaviour, we do 
not need the explicit form of p\. All we require is that the eigenvalue distribution has 
a finite extent for finite 't Hooft coupling. In particular we know that the distribution 
is non-zero only in the domain —p < x < fi. The relevant question is whether the 
solution to Eq.(3.4) is sensitive to this finite extent of the distribution. When / 

^Note that our definition of Wa^, Eq. (3.1), does not include the normafization factor given by 
the dimension of the representation. 
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is taken to be sufficiently small, we can find a saddle point with e ^'^^ ^ e^'^^, by 
expanding (3.4) in powers of e^'^^ and keeping the leading term 



^1 r dxpx{x)e-^^^ = hwu), e-2"^>e2"^ (3.14) 

J J — a J 



We already know that {Wq) ~ e^'^i^/'^ (\n \y^/'^ at strong coupling, and therefore 

the above solution applies for small / such that 



/ < e-"'^/^ ~ A'^ VMA < 1 . (3.15) 

Substituting the small / saddle point into the integral representation for the anti- 
symmetric Wilson loop, we obtain 

^ M(i^exp(7V[/ + /l„(r'(»'o))]) = (3.16) 

Note that the overall normalization given by the dimension of the representation is 
precisely cancelled by contributions to the exponent at the saddle point. Therefore, 
explicitly, we have 

^^^^^^(hUpT^' (3.17) 

Hence, the antisymmetric Wilson loop at strong coupling exhibits a crossover from 
A-dependent behaviour (3.17) for <^ A"^ vhiA to a \- independent limit (3.13) for 
parametricaly larger f = jf ■ We do not know if the cross-over is actually smooth, but 
given that the saddle-point z^p (see figure(2)) is located well away from the branch 
points, we do not expect non-analyticities in the associated configuration. Further- 
more, we have not seen evidence of more than one saddle point, so a ffist order 
transition due to competition between two or more configurations appears unlikely. 
Settling the actual nature of the crossover behaviour described above will require a 
more detailed analytical or numerical understanding of the finite-A eigenvalue distri- 
bution. 



4 Comparison with A/^ = 4 SYM 

A comparison of the result above with the corresponding one for A/" = 4 SYM is 
quite useful. Circular Wilson loops in A/" = 4 SYM are computed by the Gaussian 
matrix model with eigenvalue density 

px{x) = ^VA-(27rx)2. (4.1) 
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Using this eigenvalue distribution and upon rescaling variables in Eq.(3.6), x = ^ 
and the saddle point value of the rank k antisymmetric Wilson loop is 

determined by the condition, 



X 



2 



, di— — p = /. (4.2) 

71 J_^ 1 -|- gVHx-z) 

The crucial difference with respect to the M = 2 theory studied previously, is the 
explicit dependence on the 't Hooft coupling of the A/" = 4 theory, which now plays 
the role of the temperature of the Fermi-Dirac distribution. 

The limit of infinitely strong coupling (A — ?■ oo) can be interpreted as a "zero 
temperature" limit for free fermions and therefore all Fermi levels x below the chem- 
ical potential z are "occupied", whilst aX\ x > z remain "empty". Hence, the saddle 
point Zsp is determined by setting the total number of fermions (eigenvalues) in the 
ground state to k 

2 /"^sp k 

dyVT^ = — . (4.3) 



71 Li ^ N 



This leads to the known result [15, 17, 20], 

{WA,)M=4^^^ = exp I AT^sin^^fc I , -(^^ - sin^^ cos^^) = ^. (4.4) 



This result is only valid when A — t- oo. For any finite A, our free fermion interpretation 
above shows that there will be corrections due to occupation of all levels above 
the chemical potential. It would be interesting to look at the detailed form of the 
corrections to the infinite coupling limit. The leading corrections include powers 
of A~^ and exponentially suppressed terms ~ 6"^*-^"'^=^'^ which are suggestive of 
world-sheet instanton corrections to the D5-brane saddle point. In the free Fermi 
picture, these corrections can be thought of as contributions from excited levels at a 
finite temperature of order A"^/^. 

The salient features of the formula Eq.(4.4) are the dependences on N and y/X. 
The exponential growth with N is due to the large k limit with k ~ 0{N), while the 
factor of a/a is a consequence of the large 't Hooft couphng. In the IIB string dual 
description on AdSr, x S^, the dependence on N and a/A arises straightforwardly from 
the tension of a semiclassical D5-brane wrapping a flux supported C 5*^, which 
computes the antisymmetric Wilson loop [15-18]. The D5-brane tension T^s = 
l/{27rgs(y'^) expressed as a function of the 't Hooft coupling and AdS radius -RAds 
becomes N\/X/87r'^R\^g. Since all the curvature scales are set by -RacIS, the scaling of 
the Wilson loop is completely determined. The angular variable 6k introduced above 
determines a specific latitude in and the supergravity result precisely matches 
Eq.(4.4). 
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The result for the M = 2 SOFT suggests that the antisymmetric Wilson loop 
for large fc, in a putative large- string dual, should indeed computed by a D-brane 
since its action is 0{N) ~ g~^, where Qs is the string couphng. However, the absence 
of any dependence on the 't Hooft parameter A also indicates that this D-brane must 
probe (internal) portions of the geometry with string scale curvatures. For sufficiently 
small the antisymmetric loop is k times the fundamental Wilson loop, which, due 
to its growth with A is likely computed by a string in AdS^ with a large effective 
tension T^s = ^InA. The transition/crossover from this semiclassical picture to 
a less familiar situation could be driven by /c ~ 0{N) strings "puffing up" into a 
D-brane probing highly curved or non-geometric portions of the dual string theory. 



5 Symmetric Representation 

The symmetric tensor representations can be obtained by considering a generating 
function which is the inverse of the characteristic polynomial. 



-N / dxpxix) ln(l - e-2"(^-")) - 27ikz 



i j) dz exp 

These formulae pick out a symmetrized polynomial in the eigenvalues e^™% 

(WsAx-2 A;!(A^-1)! , r . . ^ 

^ ^ ' l<H<i2<...<ik<N 

The main difference between this and the antisymmetric tensor representation is 
that now the eigenvalues appearing in each term of the polynomial, need not all 
be distinct. In fact, we may separate out the terms in two categories: those with 
H ^ ^2 ^ ■ ■ ■ ^ ik, and those where some of the ii coincide. The contribution from 
the former is identical to the antisymmetric representation; the latter includes terms 
where most eigenvalues coincide so that their behaviour is like a multiply wound loop. 
In fact we will see below that the integral representation is dominated by two distinct 
saddle points, related precisely to the two categories of terms in the symmetrized 
polynomial. At strong coupling the saddle point related to the multiply wound loop 
grows exponentially with the 't Hooft coupling and dominates the symmetric loop. 

The integral representation above which leads to the large N limit, involves a 
function with a branch cut along the real axis —fi < z < fi. Therefore, we need 
to first determine whether putative saddle points that contribute to the integral, lie 
in the vicinity of the branch cut on the real axis. Now the large-iV saddle point 
equation can be interpreted as fixing the number of bosons at a chemical potential 
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with a density of states given by px{x), 

' dx , ^ = /■ (5.3) 

Formally, the equation for the symmetric representation can be obtained from Eq.(3.4) 
after the replacement / — > — / and z — )■ z + |. However, now the analogy with the 
Bose-Einstein distribution implies that the system could display a phase transition 
analogous to Bose-Einstein condensation as the chemical potential is dialled from 
low to high values. 

5.1 Large A and a non-analyticity 

Small / solution: We begin by locating the solution to Eq.(5.3), with the finite 
A distribution function. It is clear that the equation is likely to have solutions with 
z < — yU, for then the integrand is positive definite. A quick check will now confirm 
this expectation. Assuming that the equation is solved by some e^'^^ <^ e^^'^'^, and 
expanding Eq.(5.3) to the lowest order in e^'^^, we find a saddle point given by 



3A: 



exp(-27r2;sp) ^ (Ws,) cc -^^^^^ , / < A^^vOnA < 1 . (5.4) 

This is identical to the antisymmetric loop in the same limit, and is consistent with 
the view that in the fixed k large limit, the two should be equivalent to k coincident 
Wilson loops in the fundamental representation. In fact, we can be more precise. 
Following the same arguments that led to Eq.(3.16), we find 

i^f^\ ^ (Wo)'' for /<A-iVhrA<l. (5.5) 
dim(5fc) 

What is now interesting in the case of the sjTiimetric representation is that for a 
given (large) A, as / is increased smoothly, the saddle point moves to larger values of 
z, towards the branch-point at z = — /i. Beyond this, the nature of the saddle point 
must change as it moves through the branch cut joining z = —fi and z = +fi. Let us 
now ascertain the critical value f = fc aX which this occurs. A useful way to proceed 
is to first formally expand Eq.(3.4) as a power series in e^'^^, assuming z < — /x: 

^e'"""/ rfxpA(a;)e-'""" = /. (5.6) 

n=l ''^(^ 

The coefficients in this expansion are the VEVs of multiply wound loops in the 
fundamental representation. Following the logic of the arguments presented in [8], 
these should be determined completely by the endpoint of the eigenvalue distribution 
at strong coupling i.e. the endpoint of the Wigner semi-circle law (2.8). At strong 
coupling (/i ^ 1), we have 

Mx)e— ^i?^^, (5.7) 
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Figure 2: Left : The branch cut singularity associated to the integral representation of 
the symmetric Wilson loop lies on the real axis in the z-plane. For / less than a critical 
value fc the saddle point Zsp also lies on the real axis. Right: When f > fc the saddle 
crosses the branch point. Beyond, a qualitatively new saddle point dominates the integral: 
a single eigenvalue well-separated from the continuum, reminiscent of Bose condensation. 
A similar configuration in = 4 SYM was related by [15] to the position of a probe 
D3-brane in AdS^ x S^. 

where the right hand side follows completely from the square root behaviour near 
the endpoints of the eigenvalue distribution. This assumption yields R = 2; a. more 
careful estimate of R performed in [8], gave i? = 2.18 (for n = 1). We note that 
in principle R could depend on n, but we take this dependence to be weak. With 
increasing n, the integral is more sharply peaked near the endpoint a.t x = fi and 
we expect that taking R = 2 becomes a better approximation. The strong coupling 
saddle point equation can then be rewritten (assuming e^'^^ < e~^'^^), 

i?^Li3 (e2"(^+")) = /. (5.8) 

The critical value of / for which the saddle hits the branch point at z = —fi is, 

/, ^ i?y|c(|) A-^v^. (5.9) 

Pushing further the analogy with the boson gas, this phenomenon is reminiscent 
of the onset of Bose-Einstein condensation, as the chemical potential z is dialled 
towards the lowest energy level (eigenvalue) at x = —fi. At this point one should 
treat the lowest level (eigenvalue) and its occupation number separately, so that it 
is 'split' from the higher levels. Motivated by this line of thinking we find a result 
quite similar to that of [15] for the multiply wound Wilson loop in A/" = 4 SYM. 

Large / solution (/ > fc): Following the intuition gained from the free boson gas 
picture, we allow for the possibility that for / > fc, the occupation number associated 
to the lowest eigenvalue ai, at the edge of the distribution, is macroscopic 

This assumption alters the saddle point equations obtained from the matrix model 
(2.3) for this eigenvalue alone, while leaving unaltered the large-A^ distribution of 
the remaining — 1 eigenvalues. 
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To extract the behaviour of the symmetric Wilson loop for / > /c, we turn to 
the discrete (finite N) version of the saddle point of the matrix integral (2.3), but in 
the presence of an insertion of the generating function for symmetric loops, i.e. Eq. 
(5.1). In the large- limit (assuming rii ~ 0{N)), we arrive at two saddle point 
equations upon varying with respect to Oi and with respect to t: 



stt^ .... r . , . / 1 



Hi 



ai-K{ai)- / dxpx{x)[ K{ai-y)] = -vr— , (5.11) 

A \ai-x ) N 

The rest of the (A^ — 1) eigenvalues continue to be governed by the unperturbed 
distribution p\{x). To establish the exact form of the split eigenvalue solution for 
/ ^ /c5 just above the transition, we would need the explicit form of px{x) for finite, 
large A. However, the exact solution of the large-A^ equation (2.5) is unknown at 
finite A. Despite this, it is easy to infer the nature of this saddle point for large 
enough A or as we now see. 

We look for a solution to the above set of equations by assuming that z -C —p. 
We may then safely ignore the exponentially small contribution from the continuum 
in (5.12), so that 

^ ^ / ^ ai z + y < . (5.13) 

Using the asymptotic form of K{x) we find that the force on an eigenvalue, ai, at 
large distances from the continuum distribution, is only due to the harmonic central 
potential. This must be balanced by the constant force on the right hand side of 
Eq.(5.11), yielding. 

This result is valid as long as, 

A A; 2 , , / X 

fli ^ — /i, or — — >— InA, (5.15) 

oTT A* vr 

which is easily satisfied at strong coupling if we take / = jj ~ 0{1). Notice that 
this implies a lower bound on / for fixed large A, which is safely above fc and so the 
solution above describes the Wilson loop in the correct 'phase'. 

Substituting our solution into Eq.(5.1), and taking care to include the contribu- 
tion from the quadratic potential term in Eq.(2.3) we find 

2" 



^ exp 



2^ '^4^ 



(5.16) 



Afe 



This result has corrections to its exponent, scaling as ~ e~4jv. We have introduced 
a superscript to identify the contribution to the Wilson loop from this saddle point. 
Below we will encounter another strong coupling saddle point. 
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Interestingly, the entire analysis above could have been adapted to calculate the 
symmetric Wilson loop in A/" = 4 SYM, yielding exactly the same result as (5.16), 
provided ^ is taken to be large enough. Indeed, the known exact formula for that 
case (both from the probe D3-brane calculation [14] and the Gaussian matrix model 
[15, 19, 20]) reduces to (5.16) when \f\j^ ^ 1. These observations indicate that 
symmetric Wilson loops may be computed by a semiclassical object (D-brane) in the 
large- string dual of the M = 2 SOFT and that, in particular, the AdS^ portion of 
the geometry probed by it may not be highly curved. 

5.2 A second strong coupling saddle point 

We found above that the saddle point contribution for the circular Wilson loop in 
the symmetric tensor representation grows exponentially with A (and with A^) at 
large A. This is in stark contrast to the antisymmetric loop. However, it turns out 
that there is another saddle point which contributes to the expectation value of the 
symmetric Wilson loop and is A-independent at infinite A. The quickest way to see 
this is to note that the integrand in (5.1) and its associated saddle point equation can 
be obtained, up to an overall sign, from those for the antisymmetric representation 
after the replacements: z ^ z + ^ and / — t- — /. At this new saddle point 

W^f L^oo = (5-17) 



exp 



-N (2 In (^coshf(^,p + f) + ^) + + 27r(2sp + |)+ln4) 



coth(f(2;,p + i)) 



2 



Expanding the contribution for small f = jr 



TV' 



= exp A' 



A \ e A / 3 V A/ V A 



(5.18) 



The full symmetric Wilson loop at large A, is given by the sum of the contributions 
from the two saddle points, 

{Ws,) = W^^^ + Wf\ (5.19) 

Hence the symmetric Wilson loops can exhibit yet another type of non-analyticity, 
namely, a first order phase transition at large A due to a competition between these 
two configurations. However, for fixed / and large A, it appears that w'^'' will always 
dominate over the second saddle point, which remains exponentially suppressed. 

The existence of two saddle points for the symmetric representation can be in- 
tuitively understood from the symmetrized polynomial representation Eq.(5.2). One 
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configuration is dominated by the terms which compute the multiply wound Wilson 
loop (Tre^'^'^") and these have ii = 12 = ■ ■ - ik-, whilst the second type of configuration 
is dominated by the terms with all or most eigenvalues being distinct and is therefore 
similar to the antisymmetric loop. A similar transition has also been noted in = 4 
SYM [19, 20]. 

6 Discussion 

The matrix model of [7] provides a powerful tool for extracting predictions for a large 
class of observables in A/" = 2 SCFTs. It may also serve as a potential window into 
large- string duals of SCFTs at strong coupling. The analysis in this note presents 
us with one class of observables in the Nf = 2N theory that could shed light on the 
corresponding large- A^ dual. Below we discuss some related questions and directions 
for future study. 

Relation to Kondo models : Wilson loops in different representations can be 
viewed as impurity spins [18] coupled to an ambient theory which, in our case, is 
an A/" = 2 SCFT. These may be regarded as supersymmetric versions of the Kondo 
model (see e.g. [23, 24]). An antisymmetric Wilson loop computes the action of a 
fermionic impurity interacting with the SCFT. There is at least one SCFT namely, 
A/" = 4 SYM, for which the circular Wilson loop in the ambient theory at zero 
temperature, has been related to the impurity model (or Polyakov loop) at finite 
temperature [25, 26]. 

It is therefore not unreasonable to expect a connection between the Wilson loops 
we have computed and certain large- A^ Kondo models. In this context it would be 
interesting to understand the physical origin of the non-zero effective "temperature" 
we see in our zero temperature calculations, both for A/" = 4 theory (4.2) and for 
the A" = 2 SCFT (3.6). In the former case the effective "temperature" associated to 
the fermion impurity action scales as 1/ and parametrizes the stringy corrections 
to the D5-brane action computing the antisymmetric Wilson loop. For the M = 2 
SCFT at infinite 't Hooft coupling the corresponding equation (3.6) which fixes the 
occupation number of the impurity fermions has no free parameters and appears to be 
fixed at an effective temperature of 0{1). The resulting action for the antisymmetric 
Wilson loop Eq.(3.13) bears a striking resemblance to the impurity entropy for the 
large- A^ multichannel Kondo model of [23] (setting A' = A^ in that model). In 

''^The origin of this relation has been explained in [25] at strong coupling. In an SCFT the size of 
a circular loop, or the temperature, can always be rescaled to unity. It is not a priori clear that the 
circular loop at zero temperature should be related to the Polyakov loop at finite temperature since 
one is supersymmetric while the other has anti-periodic boundary conditions for ambient fermions 
around the thermal circle. 
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particular, both share the same logarithmic dependence on the fermion occupation 
number k/N . 

Motivated by this resemblance between the two systems, we make a purely empir- 
ical observation which appears to follow solely from the strong coupling eigenvalue 
density Poo{x) = 1/(2 cosh Suppose that we introduce a fictitious "tempera- 
ture" P into our large- equations for the antisymmetric loop so the equation for 
the fermion occupation number (3.6) reads. 



oo 



^'^''2coshf l + e2-/3(--) ^^-^^ 

In the zero temperature limit /3 — )■ oo, where only the states below the chemical 
potential z are occupied, we obtain that the saddle point (chemical potential) and 
impurity action are 

z = f Intan(^), \n{WA,)\i^^^ = -/3A^ [im Lis (i tan ^) - ^ In tan (^)] , 

(6.2) 



When expanded out for small f = jj this yields. 



' 2 iV j 36 V iV ^ ■ ■ ■ 



(6.3) 



precisely matching the impurity entropy for the large- multichannel Kondo model 
[23, 26]. The physical significance of this is unclear at the moment as we do not 
really have a means to introduce a tunable effective "temperature" in the matrix 
model for the Nf = 2N theory. In the = 4 theory, as explained earlier, the 't 
Hooft coupling appears to play this role. It is possible that studying Wilson loops 
in more general M = 2 quiver SCFT's (such as the SU{N) x SU{N) theory) might 
shed light on this connection, if any, to Kondo models. 

Aspects of large- duals: The study of supersymmetric Wilson loops in even 
larger representations, with ranks of order A^^, proved to be remarkably fruitful in 
the context of the A/" = 4 theory. The relation between the matrix model picture for 
such operators in terms of back-reacted eigenvalue distributions [15, 29], and their 
gravity duals, has been understood via explicit construction of half-BPS type IIB 
"bubbling geometries" [31, 32]. These geometries incorporate the complete back- 
reaction of the D3 and D5-branes dual to the large Wilson loop operators. It would 
be interesting to pursue the study of these large rank {0{N'^)) Wilson loops in 
the Nf = 2N theory using Pestun's matrix model. It is possible that the large- A^ 
behaviour of these operators will provide us with further clues about the large-A^ 
string dual of this theory. The first steps towards this computation are sketched in 
the Appendix. 
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Finally, it has been argued recently [10] that the string dual to the M = 2 
SCFT in the Veneziano limit is a non-critical string background with seven geometric 
dimensions. It would be extremely interesting to know if the results found for large 
Wilson loops in this paper could be naturally interpreted as expected features of 
branes in such non-critical backgrounds. 
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A 0{N'^) Young Tableaux in A/" = 2 SCQCD 

In this appendix we derive a matrix model for 'even larger' representations of SU (N), 
i.e. those whose Young tableaux are made up of 0{N'^) boxes, by restating a deriva- 
tion appearing in [30]. 

There is a well known formula for the character of an SU{N) representation, 
expressed in terms of its Young tableau: 



det(e2^(^-*)"^) 



The following formulae hold: 

det An xN = sgn{a)Y[ A a{i) (A.2) 

det(e2"(^-*)"0 = n^^'™' - ^'""' ) • 3) 

i<j 

We want to insert this into the path integral to find the matrix model for the Wilson 
loop in representation R: 

{Wr) = 1 /"[da] ■■■Tr^e^™ (A.4) 



Z 



1 



(t6Sjv « \«<i 

Manipulating the sum over the symmetric group: 

We now use the fact that the Oj's are 'dummy variables' within the integral, and so 
we can relabel them at will. This leads to 

1 r rr p2TT{N+Ri-i)ai 

- / [da] ■■■ N\ , , ^— (A.6) 

Now we label the eigenvalues. We use the notation of e.g. [30]. The eigenvalues split 
into groups labelled by the blocks of the Young tableau: 



u 



(m) 



■ ■ U 



(m.) 



U 



(m+1) 



U, 



(m+1) 
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We then have: 



/■■■^'Iin e2"(^-(^'"+2-^+*)+^^)"^" X (A.7) 



-1 

/ i<j I<J i,j 



Now we can use the identity 

/<J j,J I<J i,j 

to obtain 

/ i /<J i,j / i<j 

(A.9) 

We can again manipulate the dummy 'a' variables to write 

(A.IO) 

But using the identity for the determinant we stated earlier, the sum in the second 
term here is equal to one. Including the Vandermonde determinant, we can write 
our final result for the matrix model partition function as 

f(ii'^f'] - mii fn»"'"°-"n(«s"-4")') X (All) 

xn 



3 * 



The new terms in this expression introduced by the trace insertion are those involving 
exponentials: there is an extra force on each block of eigenvalues (labelled by J) 
proportional to Kj, and there is an inter-block interaction. 

So we write down the matrix model to be solved to evaluate large tableaux in 
the M = 2 theory: 



(Wr) 



1 — e J ' j-n (Oj — ttj \ 



X 

KJ i,j 
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This has saddle point equations: 



^af ) - - K(af )) = 1 V f - i^(af ) - a^) 1 (A.13) 



A * 2iV ^ ' ' N ^ \ jn _ Ai) 




j \ "i "i i - e » J . 

If the model behaves qualitatively like that for A/" = 4 SYM, the different groups 
of eigenvalues will be situated at intervals along the real line, with the furthest 
to the right, and each group having a profile as for the model without the trace 
insertion. This picture should be correct if one could establish that the eigenvalue 
spread within each group (~ log A) were parametrically larger than the spacings 
between groups. This would be an interesting avenue for future work. For A/" = 4 
this model yields a prediction for the IIB supergravity action evaluated on a weakly 
curved 'bubbling geometry'. The equivalent here would be informative about the 
curvatures involved in a gravity dual. 
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